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Abstract 

The KLOE detector (jl|) at DA$NE, the Frascati ^-factory, has started taking data in April 1999 
and a total integrated luminosity of 2.4 pb _1 has been collected by the end of '99, corresponding to 
~8 millions <f> decays. With these data a preliminary measurement of <f> radiative decays (f> — > 777, 
4> — > 7r°7, <f> — > 7/7, <f> — > 7r°7r°7, <p — > 7r + 7r~7, — > 777r°7 and of the hadronic decay 4> — > 7r + 7r~7r° has 
been performed. The energy spectrum of the radiated photon in case of the 7r°7r°7, 7r + 7r~7, r/-K°j 
final states allows us to extract the information on the contribution of the direct decays cf> — » /07, 
<f> — > 007. The measurement of BR( — > /07), BR((/> — > 007) can help in understanding the nature 
of /o(980) and ao(980) which is still under debate. The value of BR(<j!> — > ^'7) can be related to the 
gluonic content of the 7/(958) while the ratio R=BR(</> — » r;'7)/BR((/) — > 777) can help in establishing 
the value of the 77 — 77' mixing angle 

Furthermore a high statistics analysis of the Dalitz plot in the cj> — » 7r + 7r~7r° decay allows us to extract 
a possible contribution of the direct decay with respect to the dominant p7r mode and to obtain a 
new measurement of the parameters of the p line shape, including the p° — p mass difference. 
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1 Introduction 



1.1 Radiative decays 

1.1.1 The scalar sector: <j> — > /07, <fi — > ao7 

The lightest scalar mesons, with masses below 1 GeV, have defied classification for nearly half a century. 
The narrow states /o(980) and a (980) do not conform to standard qq quark model expectations. The 
biggest departures from theoretical predictions based on the qq model are in the total width (predicted 
r ~500 MeV, observed T ^50 MeV) and in the 77 coupling (predicted ~4.5 keV for /o, ~1.5 keV for ao, 
observed <0.6 keV for /o, ~0.2 keV for ao). Two alternative hypotheses (||) are presently under discussion 
for their nature: they could be a four quark state qqqq (R~l fm) (|J) or a KK molecule (R~1.7 fm) (|J). 
Different BR's are expected depending on their nature, as can be seen in Table [j] in the case of the /q. 



Model 


BR(0 / o7 ) 


ss( 3 Po) 


- 1CT 5 


(uu + dd)/V2( 3 P ) 


< 1CT 6 


qqqq 


~ 1(T 4 


KK molecule 


lcr 4 -^ io -5 



Table 1: Theoretical predictions for <f> — * /o7- 



Recent studies using lattice QCD (||) suggest that qqqq states occurs generically near meson-meson 
thresholds. The recent observation of 4> - * /o7> 4> ~ * a o7 with B.R. ~ 10~ 4 ([| 0) seems to be in favour 
of the qqqq scenario. The observation from Crystal Barrel of an isoscalar /o(1365) (g) and an isovector 
ao(1450) (||) as members of the 3 Po nonet allows us to search for an explanation for the /o(980), a (980) 
outside the qq model. 

1.1.2 The pseudoscalar sector: <f> — > 777, </> — > rj'j 

The reason for studying (f> radiative decays in r\ and rf is twofold since the measurement of BR(0 — > rf'y) 
can help in defining the gluonic content of the 77' while the measurement of the ratio R=BR(</> — > 7/7) / 
BR(</> — > 7/7) can help in establishing the value of the 77 — 77' mixing angle Op. 

The presence of gluon admixture in the 77' wavefunction is a longstanding problem that could be solved 
by an accurate measurement of BR(0 — > 7/7): theoretical predictions range from as low as 10 -6 in models 
with gluonium admixture or with strong QCD violations (|ll]) to 10~ 4 in different realizations of the 
quark model (jT|; H HI II)- The recent measurement of BR(</> -> 7/7) by the CMD-2© and SND@ 
collaboration at VEPP-2M seems to exclude a gluonium admixture. 

The value of Op has been discussed many times in the last thirty years ( Ell |l~9| ; p0| pi] ): the quadratic 
Gell Mann Okubo mass formula gives Op ~ —10° while experimental data give #p in the range from 
— 14°to— 20°. A recent analysis (p2|), based mainly on decays J/"0 — * VP, give #p=-16.9 ± 1.7. A crucial 
test, originally proposed by Rosner ( |l4] ) is the measurement of the ratio R=BR(0 — > ry'7)/ BR(0 — > 777). 
This ratio predicts 7.6 x 1(T 3 for P 20° and 6.2 x 1(T 3 for P 16.9°. 
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1.2 <j> -> TT+TT-TT 

About 15% of the decay into tt + tt~tt°. These final states can be due to three different mechanisms (see 
ng0): 

1. 4> — > pir decay where pn include all the three possible charge states (namely p + ir~ , p°ir° and p~ir + ) 
with the same isospin weights; 

2. <fi — ► 7r + 7r _ 7r° direct decay; 

3. e + e~ — > cj7t , with — > 7r + 7r~. 




Figure 1: Feymnan diagrams contributing to the 7r + 7r 7t final state. 



The fit of the Dalitz-plot distribution of this three-body decay allows us to discriminate between the 
three contributions. In particular it is possible to observe the direct contribution (number 2) as predicted 
by several theoretical models, that has never been observed in previous experiments. The only up to now 
published analysis of this kind ( p3] ) finds a Dalitz-plot fully dominated by the pn contribution and only 
a limit of the direct decay to be below ~ 10%. 

Furthermore a precision measurement of the p line shape parameters is also possible. In particular the 
same amount of events corresponding to the three charge states of the p, allows us to make comparisons 
between p + and p~ (CPT test) and between charged p and p°. A mass or width difference between p^ 
and p° is a signature of isospin violation as observed in other meson and baryon isospin multiplets. 

2 (J) radiative decays 

2.1 Selection criteria for radiative decays 

Some steps of the analysis and some definitions are very similar for most of the processes studied in this 
paper. 

All the processes under study are characterized by the presence of prompt photons, i.e. photons 
coming from the the LP. These photons are detected as clusters in the calorimeter [] that obey the 
relation t — r/c = 0, where t is the arrival time on the calorimeter, r is the distance of the cluster from 
the LP., and c is the speed of light. We define a photon to be "prompt" if \t — r/c\ < bat, where we use as 
time resolution of the calorimeter the parameterization <r t = 110 ps/-^/ -B(GcV). This t — r/c interval is 
often referred to as "time window" in the following sections. An acceptance angular region corresponding 
to the polar angle interval 21°-=- 159° is defined for the prompt photons, in order to exclude the blind 
region around the beam-pipe. 

1 The efficiency for photon detection is ~ 85% at 20 MeV and > 98% for energies above 50 MeV. 
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Most of the analyses described in this paper make use of a constrained fit ensuring kinematic closure 
of the events. The free parameters of the fit are: the three coordinates (x, y, z) of the impact point 
on the calorimeter, the energy, and the time of flight for each photon coming from the LP., the track 
curvature and the two angles (j> and 9 for each charged pion also coming from the LP., the two energies 
of the beams, and the three coordinates of the position of the LP. The analysis procedure adopted is the 
following: 

1. events with the appropriate number of prompt photons and charged tracks are selected from the 
"radiative stream" ; 

2. the kinematic fit is applied on these events a first time with the constraints of the total energy and 
momentum conservation and satisfying t — r/c = for each prompt photon; 

3. other selection criteria are applied to separate the signal from background; 

4. the kinematic fit is applied a second time on the surviving events with the same constraints as 
before plus other ones imposing the invariant masses of the particles present in the intermediate 
states (7r°'s, ry's etc.). 

2.2 Luminosity measurement 
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Figure 2: Comparison between data and MC of polar angle and acollinearity for Bhabha events. 





The integrated luminosity has been measured with large angle Bhabha scattering events using only 
the calorimeter information. The measurement is described in detail in ref.(^4|). The absolute error on 
the luminosity measurement can be estimated in < 3%. FigJ^ shows a comparison between data and 
Monte Carlo for polar angle and acollinearity of large angle Bhabha electrons and positrons, that are the 
relevant quantities for the evaluation of the luminosity. The agreement between the distributions is good. 
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Table 2: Analysis efficiencies for 7r°7r°7 decay and related background. 



Decay channel 


S/B 






£cut 


£tot 


Final S/B 


Signal 




76.6% 


58.4% 


88.6% 


39.6% 






0.52 


79.4% 


30.9% 


25.8% 


6.3% 


3.3 


P TT° 


3.0 


75.0% 


22.8% 


25.9% 


4.4% 


27.0 




3.1 


68.3% 


15.5% 


26.8% 


2.8% 


43.8 




0.02 


75.0% 


0.3% 


45.9% 


9 x 10" 4 


8.8 



2.3 <j> -> 7r u 7r°7 

The 7r°7r°7 final state allows us to investigate the /o7 intermediate state. The analysis scheme(^5|) has 
been developed using Monte Carlo events, simulating both the signal and the main background channels 
with the S/B ratio listed in Tab. ^. Photons coming from 7r°'s have a flat energy distribution with 
E 1 < 500 MeV while the radiative 7 is peaked at 50 MeV. The background spectrum covers the same 
energy range. The kinematic fit, applied to the 5 photon final state, is relevant in assigning the radiative 
photon and allows a partial rejection of the background. After applying a first fit with constraints on 
quadri-momentum and time of flight, the best photons' combination producing two 77 pairs with ir° 
mass and M T o T o > 700 MeV (the expected /o mass region) is selected. A second fit, imposing further 
constraints on 7r n 's mass on the assigned 77 pairs, is then performed without any assumption on f mass 
and width. This procedure correctly identifies the radiative photon in 92% of the well reconstructed /07 
events. 

After the whole fit procedure, the resulting background rejection is still not enough, especially for 
u>tt° (S/B ~ 1 after the fit). A good variable which helps in identifying the e + e~ — > ujtt process is the 
angle tp between the primary photon and the pion's flight direction in the 7r°7r° rest frame. Because of 
the different spin between /o and w ( J = 0, 1 respectively), the cost/" distribution is flat in the first case 
while is peaked at 0.5 in the second one (Fig. ^.left). Therefore the final kinematic fit is also performed 
with a different photon's assignment, requiring the best combination of four 7's into pions, which gives 
M^Oj in agreement with the uj mass. Using these criteria, the photon assignment for non-w7r° events is 
not correct and the resulting distribution is peaked at high cosi/' values for /07 and flat for the rest of 
the background (Fig. ||.left). / 7 candidates are then selected requiring cosV' > 0.8 (/o7cut) while the 
W7T° selection requires 0.4 < cos?/; < 0.8 (w7T° ut ). 

In Fig. fright, the M,OjO distribution in the /07 fit hypothesis is shown for the signal before and 
after applying /o7cut- It is remarkable that this cut is very efficient and does not significantly modify 
the M^o^o shape. In Tab. || the efficiencies for the various analysis steps as obtained from Monte Carlo 
(MC) are reported (cl: trigger, background rejection and event classification; ^ sel+x 2 : acceptance, time 
window and % 2 cuts; cut: /o7cut) together with the signal/background ratios before and after the analysis. 
Results for the W7r analysis are listed in Tab. 

The same analysis is performed on the 1.84 pb" 1 collected on December '99(p5[). Out of the 51666 

2 In the case of the /07 signal the values of the contributions to £ c i are the following: 97.9% from trigger, 90.2% from 
background rejection, 86.7% from event classification. 

3 In the case of the urn signal the values of the contributions to e c i are the following: 98.2% from trigger, 89.1% from 
background rejection, 90.7% from event classification. 
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Figure 3: Monte Carlo distributions showing the effect of the cut on cos ip Left: angle between the primary 
photon and the pion's flight direction in the 7r°7r° rest frame for ojtv (solid), /07 (dashed) and the other 
background (dot-dashed) . The fit is performed in the urn hypothesis. Right: 7r°7r° invariant mass using 
constrained variables after fitting in f j hypothesis, before (solid) and after (dashed) the /o7cut cut. 



Table 3: Analysis efficiencies for e + e — > ujtt — > 7r°7r 7 decay and related background. 



Decay channel 


S/B 






Scut 


£tot 


Final S/B 


Signal 




79.4% 


59.2% 


75.0% 


35.2% 




/o7 


1.9 


76.6% 


48.3% 


11.4% 


4.2% 


15.9 


P TT° 


6.0 


75.0% 


42.9% 


40.5% 


13.0% 


16.2 




6.0 


68.3% 


17.2% 


23.9% 


2.8% 


75.4 




0.04 


75.0% 


0.5% 


43.1% 


1.7 x 10~ 3 


8.3 
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Figure 4: cosip variable before applying W7r° ut (left) and ir°j invariant mass after this cut (right). Fit 
constraint results have been used. In green the expected background contributions from Monte Carlo 
is superimposed while solid histograms are simulated distribution including signal summed up with the 
background. 



events with at least 5 neutral clusters, 37678 are in the acceptance angular region and are reduced to 
1815 after applying the time window requirement. Performing the u>ir° fit, 980 events have a good x 2 - 
The distribution of the cosip variable for this sample is shown in Fig. [|.left together with the expected 
one from Monte Carlo. Because of the excellent agreement, the ^7r[? ut is performed. A nice peak at the u> 
mass appears in the M 7r o 7 distribution of the surviving 529 events (Fig. fright). Background evaluation 
from Monte Carlo gives a final signal counting of 436 ± 25 (stat). Assuming a global systematic error of 
8%0 and correcting for luminosity and analysis efficiency we quote a cross section 



cr(e + e 



7r u 7r u 7) = ( 0.67 ± 0.04 (stat.) ± 0.05 (syst.) ) nb 



(1) 



in good agreement with the SND measurement (| 

The /q7 fit applied to data yields 679 events with a good x 2 out of which 307 survive also the /o7cut- 
In Figs, ^.left (center) the radiative photon's energy and the invariant mass of the 7r°7r system are shown 
after kinematic fit for the events before (after) applying this cut. In the same distributions the expected 
background contribution, estimated from MonteCarlo, is reported in the solid coloured shapes; the S/B 
ratio after the cut improves of at least a factor 3, as expected, and a clear peak above background appears 
around 950 MeV in the invariant mass (right). 

After subtracting the 112 ± 11 background events, the total counting for the signal is 195 ± 20 (stat.). 

4 The preliminary estimate of the systematic error has the following contributions: 5% from classification, 2.5% from 
clustering (50% of the events with at least 1 cluster not correctly reconstructed), 4% from fit (50% of the wrong MC 5 
photon's assignment) and 3% from luminosity. 
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Figure 5: Distributions of the radiative photon's energy and of the invariant mass before (left) and 
after (center) the /o7cut- The background contribution, estimated from Monte Carlo, is superimposed. 
The green region represents ojtt events while in yellow all other contributions are summed up. Right : 
same distributions after background subtraction. 



As in the urn case, 8% systematic error is assigned to the measurement. We obtain for M^n > 700 MeV: 



BR(0 -> /q7 vr 7r°7) = ( 0.81 ± 0.09 (stat.) ± 0.06 (sist.) ) x 10~ 4 



(2) 



to be compared with the results of the Novosibirsk experiments^; [?]). 

This measurement is still to be considered preliminary since work is in progress to estimate the 
systematic errors directly from the data and to correct the analysis efficiency as a function of in 
order to calculate a BR independently from any Monte Carlo assumptions. 



2.4 (j) -> /o7, with / -> tt+tt 

The analysis of the <fr —> /o7 decay in the charged channel /o — > tt + tt~ has been performed on a sample 
of 1.8 pb -1 of collected data by looking at the spectrum of the production cross section of ititj events as 
a function of 7r + 7r~ invariant mass squared, Q 2 (p7j). Two other processes contribute to the 7r + 7r _ 7 final 
state: Initial State Radiation (ISR) , in which the photon is emitted by the incoming electron or positron, 
and Final State Radiation (FSR), in which the 7 is emitted by one of the two pions. The latter process 
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gives rise to an interference with the signal whose sign is not known. 

The 7r + 7r _ 7 events are selected using both drift chamber and electromagnetic calorimeter informa- 
tions. The first step of the signal selection requires a prompt neutral cluster and a vertex close to the 
interaction point. This general selection identifies not only 7r + 7r~7 events, but also /x + fi~j events and a 
huge amount of radiative Bhabhas. 

The kinematical properties alone are not enough to suppress the eej events, therefore a likelihood method 
has been developed which uses both the particle time of flight and some informations coming from the 
cluster associated to the particle. This method has a 95% selection efficiency for a pion and a ~ 94% 
rejection power for electrons. 

After the likelihood selection, kinematical cuts have been chosen in order to get a further reduction 
of /i + /i~7 and e + e~7 background and to emphasize the 4> — ► 7r + 7r~7 decay contribution. For both these 
purposes particles have been selected in the central part of the detector (45° < 6 < 135°), since the polar 
angle distribution of the charged tracks from e + e~j and /i + /i _ 7 events and of the photon from ISR are 
enhanced at small angles. 

The last cut selects events based on the invariant mass of the charged track identified by applying 
4-momentum conservation in the hypothesis of a massless neutral particle: 

ifl + P2? - (m$ - y/ft + Atf R - y/% + MQj = 

where p± and p2 are the tracks momenta and Mtr is the track mass, assumed to be the same for both 
charged particles. The distribution of the variable Mtr before and after the likelihood cut is shown in 
figj|. The pion peak is clearly visible and the signal events are selected in a window of ±10 MeV around 
M n = 140 MeV, the central value of the fit. 

In order to compare the experimental results with the theoretical predictions for e + e~ — * 7r + 7r~7 pro- 
cess, the Q 2 spectrum has been corrected by the total selection efficiency as a function of Q 2 , represented 
in fig. [?]. The low efficiency at high Q 2 values is due to the trigger cosmic veto, which can mistakes a 
7r + 7r~7 event with a soft photon for a cosmic event. 

At low Q 2 values the efficiency decreases because of the cuts applied to reduce tt + tt~tt° background. 

Owing to the high rate of <p — > 7r + 7r~7r° decays, 7r + 7r~7r° background can survive 7r + 7r~7 selection: the 
contamination is higher at low Q 2 values, and it has been evaluated by fitting the track mass distribution 
of events identified as 7r + 7r~7r° for each Q 2 bin and extrapolating the fit function in the ir + ir~j mass 
region. 

The experimental data, corrected by the total efficiency, have been fitted using the theoretical spec- 
trum including ISR and FSR contributions only. The fit parameters are the integrated luminosity and the 
normalization factor of 7r + 7r~7r° spectrum with respect to 7r + 7r~7 one. In fact to evaluate the ir + ir~ir° 
background contribution to the experimental differential cross-section, for 7r + 7r~7r° events has been as- 
sumed the same selection efficiency as for 7r + 7r~7 ones, the latter fit parameter (Norm) takes into account 
possible differences in the overall efficiency. 
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Figure 6: Distribution of the variable Mtr (see text) before and after the likelihood selection. The pion 
and muon peaks are clearly visible. 
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Figure 7: Total selection efficiency for a 7r + 7r~7 event in the central part of the detector; the trigger 
cosmic veto determines the low efficiency at high Q 2 values, while the low Q 2 inefficiency is due to the 
kincmatical cuts applied to reduce 7r + 7r~7r background. 
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Figure 8: Experimental cross-section as a function of Q 2 compared to the theoretical one for pure QED 
contributions. 



The results of the fit are: L = 1774 ntT 1 , Norm = 1.16, with a x 2 /ndf = 33/45. Fig.| shows 
the comparison of the experimental differential cross-section for the e + e~ — > 7r + 7r _ 7 process with the 
theoretical one. The comparison is good and no fo signal is needed, with the available statistics, to 
describe the spectrum. An upper limit on the branching ratio for the decay <f> — > /07 — > 7r + 7r~7 can 
be set by fitting the Q 2 spectrum in the region Q 2 < 0.84 GeV 2 and extrapolating it in the photon 
energy range 20 MeV < E 1 < 120 MeV, where the signal is expected. An excess of 35 ± 160 events, with 
respect to the ones predicted by pure QED, is found. Assuming the isospin symmetry and ignoring the 
interference with FSR, this number corresponds to an upper limit on the value of the branching ratio of: 

BR(0 -> /07 -> 7r+7r- 7 ) < 1.64 x 10~ 4 @ 90% C.L. . 
2.5 <p — > ^7r°7 with Tj — > 7T 

This process is characterized by 5 prompt photons without charged tracks in the final state. It is expected 
to be dominated by the <fi — > ao7 decay, in which the ao(980) decays into rjir . The spectrum of the photon 
radiated by the 4> is expected to be broad and peaked at ~ 50 MeV. Two other processes contribute to 
this final state: <f> — > p°7r° and the non-resonant process e + e~ — ► lot: with the rare decays of p° and uj 
into 777. 

The main background comes from 7r°7r°7 final state, which is dominant in the 5 photon sample; the 
expected number of events is 10 times bigger than the signal. The other relevant background comes 
from the <f> — > 777 decay, with 3 and 7 photons in the final state, that can be reconstructed as 5 photon 
events due to accidentals in the calorimeter or photon splittings and mergings. According to the MC the 
probability for both processes to be reconstructed as 5 photon events is about 3%, then due to their high 
branching ratio the expected number of events in the 5 photon sample is also of the order of 10 times the 
signal. 
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The events are selected by requiring: 

1. no tracks in the drift chamber, 

2. total energy in the calorimeter greater than 900 MeV, 

3. exactly 5 prompt photons in the angular acceptance region. 

On the 2200 events selected in the 2.4 pb _1 sample, a first kinematic fit has been applied by imposing 
constraints on the total energy and momentum conservation, and on the consistency of time and position 
in the calorimeter (t — r/c = 0) for each photon. A cut corresponding to P(x 2 ) < 1% has been applied. 

Then for each event three different variables are constructed in order to test the three hypotheses: 



1. rpfl-y hypothesis: D m o, = liMl^M^l ^ iM^-M^ 



o -L • n (Mi2-«,o) ! . (M34-AL0) 2 

2. 7r u 7r u 7 hypothesis: £' 7r o 7r o 7 = . 3 g + - 3 J 



>r 



3. m hypothesis: D„ = J «E + ISEgffi 



The value of each £)- variable is obtained by choosing the photon pairing that minimizes it. M12 and 
M34 are the invariant masses of the photon pairs, E ra d=363 MeV. <r v — 20 MeV and er^o — 9 MeV 
have been evaluated from the data themselves, by fitting the invariant masses distribution on a sample 
of events. u ra d is obtained from the energy resolution of the calorimeter. 

The following cuts are applied: 

1. D„ n o~ < Z3 7T o 7r o 7 in order to select the events that have a bigger probability to be r)ir 0r y rather than 
7r°7r°7 (see fig.^) 

2. D vl > 2; in fig.|l(] is shown that this cut is able to reject events in the peak of E ra d at 363 MeV, 
that correspond to the 77 mass peak. 

On the 240 events selected a second kinematic fit is applied, imposing as further constraints the two 
invariant masses of r\ and ir°; a cut corresponding to P(x 2 ) < 1% h a s been applied, and the 74 surviving 
events form the final sample. 

The efficiencies for signal and backgrounds have been evaluated by MC, taking into account the 
dependence on photon energy. The resulting efficiencies are listed in Tab.^: in the first column is reported 
the trigger plus the background filter one; in the second column, the efficiency of the 5 prompt photon 
cut, that is mostly due to the angular cut at 21°, and in third one is reported the effect of the selection 
and of the two kinematic fits. 

In fig.|ll] are reported the invariant mass spectrum of the four photons assigned to r\ and 7r0 and 
the distribution of the cosine of the polar angle of the radiated photon, which agrees with the expected 
1 + cos 2 9. 

The background spectrum, superimposed in fig.|ll|, is a prediction obtained by MC by weighting the 
background processes with the product of their expected branching ratio times their efficiency. 

In order to evaluate the Br(<f> — > i]ir 0r y) we consider the whole spectrum: 74± 9 events are selected 
(N), with an expected background of 21± 6 (B), and assuming Br(rj — > 77) = 39.2% 

Br(<p -> r/7r° 7 ) - ^ ~ = (0.77 ± 0.15 stat ± 0.11 SJ/st ) • 10" 4 (3) 

eLa^Brir) -> 77) 
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Figure 9: £>-variables for the r)ir°~f and 7r°7r°7 hypotheses; a cut is applied on their difference by selecting 
the positive part of the distribution. 




Figure 10: Rejection of 777 background, the circle in the scatter plot shows the rejected region with the 
cut D m > 2. 
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Process 


Trigger + 


5 prompt 


Selection + 


Total 




niters 


photons 


kin. fits 




<j) — >• a 7 -> ?77r 7 


0.81 


0.70 


0.40 


0.23 


— > /9°7T — > r77T°7 


0.77 


0.70 


0.26 


0.14 


<j>-> fol -> 7r°7r°7 


0.80 


0.70 


2-10" 3 


10~ 3 


p°7r° -» ^°7r°7 


0.80 


0.70 


0.03 


0.02 


e + e~ — > W7T° — > 7r°7r°7 


0.82 


0.70 


0.02 


0.01 


<j> -> r]>y 


0.78 


0.03 




< 5 • 10~ 4 



Table 4: Efficiencies evaluated by means of the MC simulation; for the e + e — > uitt° — > 7]7r°7 process we 
assume the same efficiencies of the 7r°7r°7 final state. 
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Figure 11: Upper plot: invariant mass of the rjir system, histogram: data; shadowed histogram: expected 
background (from MC); black points: data - background difference. Lower plot: distribution of the cosine 
of the polar angle of the radiated photon, black points: data; histogram: MC - signal only. 
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where L = 2.4 pb _1 , and = 3.2 fib. 

This value is in agreement within the errors with the results of the Novosibirsk experiments^; |(]). 

According to the MC we expect 6 events from the processes <f> — ► pn — > rjir°j and e + e~ — * ujtt q — ► 
777r°7, then considering them as background we obtain: 

£r(</> a o7 -> ^° 7 ) = (0.69 ± 0.14 stat ± 0.10 sys4 ) • 10~ 4 (4) 
2.6 — > t/7, with T) — > 77 

The </> — ► 777 — * 777 decay, having an higher BR (0.49%) and an harder energy spectrum with respect to 
the rest of </> radiative decays, is a good calibration sample for multi-photon final states. 

Since the energy of the radiative photon (E lrad ~ 360 MeV) is inside the energy range of the 7's 
coming from the 77 (150 < E 7n < 500 MeV), the kinematic fit is used to select the 77 pair assigned to the 
meson. For each event with three photons in time window the fit is applied in the 777 hypothesis using 
the 77 mass constraint for the three possible photons' combination. The minimum \ 2 is then selected. 

The only relevant background for this channel comes from the <j> — > ir°j decay (S/B ~ 4). Since 
the energy of the primary photon (E liad ~ 500 MeV) is higher than the one of the signal, the fit can 
reconstruct it as an 777 event by wrongly combining the radiative photon with a 7 coming from the 7r° 
(E 1 < 500 MeV). Being the third photon constrained in the 360 MeV region, the three photons are 
monochromatic for 7r°7 events after the fit. The background is therefore rejected using the AE variable, 
energy difference of the 77 pair, requiring |A.E| < 330 MeV: as it is shown in Fig. |l2|. up-left, 7r°7 events 
are peaked at high \AE\ values while signal has a flat distribution. 

A sample corresponding to 1.84 pb _1 have been analysed: among 226736 events with at least 3 
clusters, 183345 satisfy angular acceptance cut and 45889 have also 3 photons in time window with 
E > 20 MeV. After the \ 2 cu t an d the background rejection the sample is reduced to 18504 events. The 
resulting angular and energy distributions of the radiative photon, together with the 77 pair invariant 
mass, are shown in Fig. |l2|. All distributions are well in agreement with Monte Carlo expectations. The 
resulting 77 mass is within 0.2% the expected value. 

In order to check if 777 QED background can simulate the signal, angular distributions between 
photons' pairs have been studied (Fig. |l3| ). The excellent agreement with 777 simulated events leaves no 
room for any residual background. 

The <j) visible cross section has been evaluated by measuring 

a{e+e- ^0^777^ 777) = ^ 7 • — — r (5) 

L x e ana BR(0 -> 777 -> 777) 

and using BR(0 — > 777 — > 777) = (0.49 ± 0.02)%. Luminosity is estimated by using large angle Bhabha's 
(6 > 45°) while analysis efficiencies, listed in Tab. || ^] are evaluated from Monte Carlo. 
We quote: 

o (p = ( 3.19 ± 0.02 (stat.) ± 0.23 (syst.) ) fib (6) 

This value is in good agreement with CMD-2 result, obtained using K$ — > 7r + 7r~ events: — (3.114 ± 
0.034 ±0.048) Mb®. 

Complete evaluation of systematics is in progress. At the moment a preliminary estimate of the 
contributions is summarized in Tab. |[ 

5 The values of the contributions to £ c j are the following: 92.6% from trigger, 81.5% from background rejection, 97.0% 
from event classification. 
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Figure 12: Comparison between data ( — ) and Monte Carlo (•) distributions for <j> — > 777 — > 777 events: 
energy difference of the two photons assigned to 77 (up-left) - the two peaks are due to the cj) — > ir°j 
background; angular distribution (up-right) and energy spectrum (down-left) of the radiative photon; 77 
pair invariant mass (down-right). Last two variables are obtained after photons' assignment from the fit 
but using energies reconstructed without constraints. 
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Figure 13: Opening angles between photons' pairs of — > 777 — > 777 events for data and Monte Carlo. 
Photons named 1 and 2 are the ones assigned to r\. 
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Table 5: Analysis efficiencies for <j> — ^ 777 — > 777 events. 



Efficiency 


Contribution 


e cl = 76.0% 
£ sc i = 93.7% 
e x 2 = 98.5% 
SAE = 91.3% 


Trigger, background rejection and event classification 
Acceptance and time window cuts 
X 2 cut 
AE cut 


e ana = 64.3% 


Total 



Table 6: Error contributions to cross section measurement using (f> —* 777 — > 777 events. 



VARIABLE 


VALUE 


METHOD 




4% 


From PDG 


L 


3% 


Comparing trigger vs L3 vs offline values 




5% 


Conservative estimate using a small data 
control sample without event classification 


£scl 


1.0% 


Fit to cos# 7rac i distribution in different angular 
regions for acceptance. No contributions from 
TW cut (TWmax 2 — > 3 ns: no changes) and 
clustering (emc — 99.9%) 


e x 2 


0.6% 


50% of the wrong assignments percentage as 
obtained from MC 


SAE 


0.3% 


Comparison between loss of events counted 
with MC and the estimate on data fitting the 
\AE\ < 300 MeV region 
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In addition to the analysis described above an alternative method has been developed for measuring 
the ratio Br(0 — * j]j)/Hi((j> — > ir°j) without using any kinematic fit([29"|). 

2.7 The ratio BR(0 -> 7/7 ) /BR(</> -> 777 ) 

Two different decay chains giving rise to both fully neutral and charged/neutral final states have been 
used to study the ratio R$ — BR{<j) — ► 7/7 )/BR((j> — ► 777 ): 

1. ^ — * 7/7 — * 777r + 7T~7 — * 7T + 7r _ 777 

2. <^> — * 7/7 — + 7777° 7t° 7 — * 77 

A very clean control sample is given by <fi — > 777 decays with identical final state: 

1. (f> — ► ?77 — * 7T + 7r~7r° 7 — ► 7T + 7T _ 777 

2. — > 777 — > 7T° 7T° 7T° 7 — > 77 

These events, being 2-3 order of magnitude more probable than the corresponding <j) — > 7/7 ones, consti- 
tute also the main source of background for their detection. For this reason a kinematic fit with mass 
constraints is needed to obtain a satisfactory signal to background ratio in these final states. 

Since the final state is identical for the <p ~ * v'l an d 4> ~ > VI corresponding channel, most of the 
systematics will cancel if we evaluate the ratio using the same final state to count <j> — > 777 and 
4> ~~ * v'l events. 

2.7.1 7r + 7r _ 777final state 

For <fi — * rf'y events this final state is characterized by a nearly monochromatic photon with E 1 = 60 MeV 
recoiling against the 77', and two (harder) photons coming from 77 — > 77 annihilation. On the contrary, for 
<f> — > 777 events the radiative photon, still monochromatic, is the most energetic one, with E 1 — 363 MeV 
and the two (softer) other photons come from 7r° — > 77 annihilation. 

In addition to the <f> — > 777 background, some background events can be expected from cj> — » K$ Kl events 
with one charged and one neutral vertex where at least one photon is lost and the Kl is decaying near 
the interaction point and from <f> — > 7r + 7r~7r° events with an additional cluster counted. 

A first level topological selection runs as follows: 

• 3 and only 3 prompt neutral clusters (as described above) with E 1 > 10 MeV and 21° < # 7 < 159° 

• 1 charged vertex inside the cylindrical region r < 4cm; \z\ < 8cm 

and is common for both cj> — > rf'y and cj> — » 777 events. 

Background from 7r + 7r~7r° events is strongly reduced by means of a cut on the sum of the energies of 
the charged tracks assumed to be pions: it is expected larger for three pions events than for radiative 
events. 

Background from (f> — > K$ Kl is reduced for — > 777 using the fact that the spectrum for photons 
coming from kaons is limited to energies below 280 MeV while in cp — > 777 we expect at least one photon 
with energy exceeding 300 MeV. 

The same cut cannot be applied to (f> — > r]'f events where the energy spectrum of photons is different: 
in this case however, a suitable variable to select the signal has been found to be the sum of the energies 
of the three photons E~ lll . In conclusion one applies the cuts: 



19 



• For (f> — > 7/7 selection: 

- 7^+ + E v - < 550 MeV (e^ ~ 1.5 • 10~ 3 ) 

- ^ nax > 300 MeV (e^ Ks - 2 • 10~ 4 ) 

• For (p — > r]'"f selection: 

- E„+ + E^ < 412 MeV (e 37r ~ 1 • 10~ 4 ) 

- £ 777 > 520 MeV {sk l Ks - 1 ' 10~ 4 ) 

The efficiency for this selection is evaluated from Monte Carlo to be 39.6% for tj> — > j/7 and 37.9% for 
— > 777 events. 

Contamination from — ► 777 events into the <^> — * 7/7 sample is at this level still very high (S/B 
~ 10 -3 ): thus a kinematic fit with mass constraints (see section 2.1) has been implemented for both the 
decay chain hypotheses constraining all intermediate masses. The energy spectrum of the photons gives 
no problem in assigning clusters to particle originating them: as already noticed above radiative photon 
is the most energetic one in <fr — > 7/7 events, while it is the less energetic one in <fi — > 7/7 events; the 
other two cluster belong to n° and 77 respectively in the two cases. A cross cut on V(xi, ) > 25% and 
V(Xm) < 1% has proven by Monte Carlo to maximize the significance S/yfS and has thus been chosen 
as final selection criterium for 4> — * 77' 7 events. Final Monte Carlo efficiency after this cut is 18.6% for 
4> — > rj lr y while for 4> — ► 777 a 90% C.L. upper limit can be set to 4.4 • 10 -5 giving rise to an expected S/B 
ratio > 35 (90% C.L.) if one uses the PDG'98 value for BR(0 -> 7/7 ). 

A selection cut can also be put to select <fi — > 777 events, and has been chosen in a very conservative 
way to be P(Xn 7 ) > 1% due to the low background on this channel. With this cut one has a final 
efficiency of 31.9% for cf> — > 777 and selects a very pure set of events with background (estimated from 
Monte Carlo, and confirmed by a fit to the 77 mass peak) being ~ 0.1% of the sample. 

The abundant and pure <f> — > 777 events can be used as control sample for systematic effects on the 
efficiency, and to compare data versus Monte Carlo distributions for the variable on which the cuts are 
set. 

All comparisons (see fig.|l4]) show very good agreement (within 1-2%) between data and Monte Carlo, 
and since the dependence of the efficiency on the cuts is not critical (for example moving the cut on 
the charged pions energy by ±1% changes — > 777 selection efficiency by ~ 0.1%, for a more detailed 
discussion see (p0|)) the overall systematic error on the estimation of efficiencies is very small. Also, 
when evaluating the ratio most of the systematics will cancel out due to the strong similarities 
between the two categories of events. With the statistics of ~ 2.4 pb -1 of 1999 run we found 21± 4.6 
4> — * ffl events in this decay chain with less than one event of background expected at 90% C.L., while 
with the <j) — > 777 selection selects 6696 events in the same runs. The distribution of the invariant mass of 
the two charged pions and the two most energetic photons in the event is shown in fig. ^| compared to 
the Monte Carlo expected for pure cp — > 7/7 events. Solving for R$ we get: 

N v , y e vl BR(ri -> tt+tt-tt )BR(tt° -> 77) 



R 

and, thus: 



V„ 7 e V7 BR(ri> -> 7^-77)7^(77 -> 77 ) 



7? = (7.1 ± 1.6(stat.) ± 0.3(syst.)) ■ 10~ 3 

where the systematic error is dominated by the uncertainty on the value of the intermediate branching 
ratios (4%). In fact systematic effects on luminosity evaluation and cancel out exactly in the ratio, 



20 




Figure 14: Monte Carlo (pure <j> — > 777 ) versus data for tt + tt 777 events: a) Cluster energy spectrum; 
b) Charged pions momentum spectrum; c) r\ invariant mass distribution; d) x 2 of kinematic fit 
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Figure 15: Invariant mass of two charged tracks and two most energetic photons for events selected as 
4> — * v'l ■ Data (lower plot) is compared to pure monte Carlo <j> — > r]'j events. 
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while trigger efficiency, streaming efficiency and reconstruction effects cancel almost exactly. This result 
on i?0 leads in turn to: 

BR{<f) -> 7/7 ) = (8.9 ± 2(stat.) ± 0.6(syst.)) • lCT 5 

This result has to be compared to the most recent results by CMD-2 and SND, which use the same final 
state (H|; |§) 

2.7.2 77 final state 

The high hermeticity of the KLOE EmC allows us to detect with high efficiency multiphoton final states. 
The 4> ~ > vfl ~ * 71-0 7r ° 7 77 — > ?7 anf i its background <p — > 777 — > ir° ir° ir° 7 — > 77 are two such states, and 
can be selected by seeking seven prompt clusters in the EmC with no charged track in the event. 
As far as other backgrounds are concerned, only fully neutral channels are relevant. However, due to 
high hermeticity of the EmC, the most relevant of them, the cf> — > Kg Kl — > 5ir° — ► IO7 , may mimic 
a 77 final state only in a very small fraction of events. Moreover, the prompt 7 selection rules out all 
events where Kl is not decaying very near the beam pipe. A sample of 5 • 10 5 Monte Carlo events has 
been generated for this background, and no event survived the topological selection cuts. The simple, 
topological selection is common to both <f> — > 777 and <p —> r]'j events and requires: 

• 7 prompt neutral clusters with 21° < 9 < 159°. 

• No charged tracks. 

• \E tot - 1020 MeV| < 130 MeV 

where E to t is the sum of the energies of the selected clusters. 

The efficiency of this cut is 41.3% for <f> — > 777 events and 41.2% for <\> — » 7/7 events. Events passing this 
cut are further analyzed in both hypotheses of being <p — > 777 and <f> —* rj'j events. 

To completely rule out any cf) — > K$ Kl background a cut on E™ lax > 300 MeV can be put for 
4> ~ y Vl events with essentially no loss in efficiency. Monte Carlo simulations, and fit to the 77 mass peak, 
show that background in the <j> — > 777 channel is then again at the level of 0.1%. 

No attempt is made to solve the combinatorial for the three ir° 's coming from the 77 decay, and the 
only identified photon is the radiative one, which is, as usual in cf> — > 777 , the most energetic photon of 
the event. 

For <p — > 7/7 the further analysis is based on a two-steps kinematic fit. First, a kinematic with no 
mass constraints is performed to achieve a better determination of the photons energies. Then a pairing 
procedure is applied in order to obtain the correct identification of the photons. 

The starting point of this procedure is the observation that the most energetic photon of the event 
comes always from the 77 decay. The remaining six photons are then scanned to check the best pairing 
giving the correct 77 mass. Once the second photon is assigned to the 77, the most energetic of the five 
remainders is coming from one of the tt° 's : Monte Carlo shows that in this way the first three photons 
are correctly assigned in 99.3% of the cases. 

The remaining 4 photons give rise to twelve possible combinations: a % 2 -like function is built for 
each combination to compare the obtained masses to the expected ones, and among these the best five 
combinations are selected : in 96% of cases among these there is the correct one. Finally the five best 
combinations are fitted with a full constrained kinematic fit, where the 7r° 's , rj and 77'invariant masses 
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Figure 16: Monte Carlo energy spectrum of the most energetic cluster for events passing both <j) — > 
777 and <f> — > j/7 event selection, and after the kinematic fit with no mass constraints. The plot shows 
4 1 ~ > 11 events (lower plot), <p ~ * I'l events (upper plot) and the position of the selection cut (arrow). 



constraints are added to the constraints used in the preliminary fit. The combination minimizing the 
final x 2 is then chosen to be the correct one, and the x 2 °f this fit is used as a discriminating variable 
for background suppression. 

The combination found is the correct one in 90% of the cases, the main sources of mistakes being the 
radiative 7 associated incorrectly to a tt° (5%) and/or the 4 7 's of the tt° 's being mismatched (5%). 

The x 2 corresponding to the best combination is compared to the one of a kinematic fit performed 
in the 4> — > 777 hypothesis, with only the r\ mass as intermediate mass constraint. Using Monte Carlo 
events an optimized cut has been chosen in a triangular shaped region in the plane V(x 2 >~ ( ) — 'P{x 2 -y)- 
Analytically it can be described by the formula: 

Vix 2 ^) > 15% + 3.4-7>(x" 7 ) 

This cut alone, although being able to reduce drastically the 4> —> VI background, is not enough to obtain 
a satisfactory Signal/Background ratio: indeed 5 x 1CP 4 <f> — » 777 events still survive the cut against 15.6% 
of <f) — > v'l ' giving a S/B ratio of ~ 0.8. 
For this reason a further selection cut has been introduced. 

The distribution of the energy of the most energetic photon after the kinematic fit is shown in fig. ^| 
for <j) — > 7/7 and <f> — > 777 fully neutral events. A cut on this energy is able to scale down definitively 
the ^ — > 777 background, even if it causes a somewhat loss in efficiency for the — ► rj'j signal. The 
maximization of the significance lead to the choice of a cut at E^ ax < 340 MeV. 

The final selection efficiency for <f> — > 7/7 events is £s e i = 12.7% while selection efficiency for the 4> ~ ¥ 
777 background goes to ~ 6.7 x 10~ 6 . This results in an expected S/B ratio > 20 (90% C.L.) using the 
PDG'98 value for BR(0 -> 7/7 ). 
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Figure 17: Monte Carlo (pure cj> —> VI sample) versus data for 77 events: Cluster energy spectrum (left); 
•q invariant mass (right) 

Applying the selection criteria described above to the 2.4 pb _1 statistics of 1999 runs, we select 
6^2 2 ^ ~ * v'l ~ * ?7 events (with less than one event of background expected at 90% C.L.) and 10938 
4> ~^ "HI ?7 events. This is the first observation of the decay chain <fi — > 7/7 — > 77. The distributions 
of the <f> — > 777 control sample compare very favourably with the simulations for what the variables on 
which the cut are set are concerned (see plot |l7|). This gives, for i?^: 

_ N, Yl Srn BR (V 71-0 7r ° 7r ° ) BR ( 7r ° zl 77 ^ 
~~ N V1 e v >~, BR(rj' -> 7T° tt° ri)BR(v ~> 77 ) 

and, thus: 

R 4 = (6.911:5 (stat.) ± 0.9(syst.)) ■ 10~ 3 

where, as already discussed for the 7r + 7r - 7 77fmal state, since most of the systematic effects cancel 
out in the ratio, the systematic error is dominated by the efficiency for low energy photons (10%), by 
the uncertainty on the value of the intermediate branching ratios (6%) and by a 5% systematic effect 
evaluated on the efficiency of the x 2 cut due to radiative photon misassignment. It is interesting to note 
that the analysis performed in this channel, although being statistically less accurate than the one for 
charged/neutral final state, is fully compatible with that one, and constitutes the first measurement of 
Re], performed with a decay chain different from the one leading to the 7r + 7r~7 77final state. 

2.7.3 The 77 — 77' mixing angle 

The importance of the measurement of R^ to extract with great precision the pseudoscalar mixing angle 
has been stressed many times during the years (|l4|). If one neglects <fi — uj mixing and SU(3) breaking 
effects in the effective Lagrangian, the <fr — > 777 , 4> ~ v v'l decays may be described as simple magnetic 
dipole transition, giving for the ratio R^ the value: 

= cot 2 <p P (^j 

with ipp = 'dp + arctan\/2- In a recent paper by Bramon et al.(p2|) it has been stressed, however, that if 
one takes into account also <fi — oj mixing angle tpy — +3.4°, and accounts for SU(3) breaking via a term 
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Parameter 


Fit result 


PDG result 


M( Po ) (MeV) 


776.1 ± 1.0 


776.0 ±0.9* 


AM (MeV) 


-0.5 ±0.7 


0.1 ±0.9 


r(p) (MeV) 


145.6 ±2.2 


150.9 ±2.0 


A(direct term)/A(p7r) 


0.10 ±0.01 


-0.15 -=-0.11 


fase(direct term)-fase(p7r) 


(114 ± 12)° 





Table 7: Results of the fit to the Dalitz plot compared to the PDG values, 
proportional to ^ ~ 1.45 the formula above gets a correction factor: 

\ m sm Zipp J \p v J 

The formula above has been used, together with the measured value of to extract a measurement for 
■dp , giving: 

$p = (-18.9 o ±f;go(stat.)±0.6°(syst.)) 

3 (j) — > 7r + Tr _ 7r° 

7r + 7T _ 7r° events are selected requiring a prompt vertex with two opposite sign tracks and two prompt 
photons in the calorimeter. From the two tracks, the direction of the missing momentum can be evaluated 
and associated with the ir° direction. The opening angle between the two photons in the ir° rest frame 
is required to be larger than 170°. Furthermore in order to remove a residual background mainly due to 
e + e~77 an opening angle between the two tracks less than 170° is also required. 

The final sample (330000 events) has been analyzed by means of the Dalitz-plot method. 

The Dalitz-plot binned in 8 x 8 MeV squares and corrected for the efficiency is fitted to a model of 
7T + 7r~7r° production including the following terms: 

1. A P1T is the pir amplitude given by the sum of the three p charged states, each described by a 
Gounaris-Sakurai parametrization. Free parameters are the p masses and the width; 

2. Adirect is the direct term contribution given by a complex number that is two free parameters, 
namely a modulus and a phase. The modulus is normalized in such a way that a value equal to 1 
corresponds to a direct term equal to the pir term. 

3. A^Tr is the lot: term, where mass and width of the tu are fixed to the PDG values, and only a complex 
amplitude that is a modulus and a phase is let free. 

The fitting function is then given by (X and Y are two Dalitz variables): 

f(X, Y) = \p + x p -\ 2 ■ \A p7r + A direct + A^| 2 

where the square of the vector product in front takes into account the vector nature of the decaying 
particle. 

In Table.0 the results of the fit are shown and compared with PDG values. Two observations can be 
done. 
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First we observe a sizeable direct term (about 10 % of the pw term) with a phase respect to pit loosely 
close to 90°. We remark that this is the first observation of this decay. 

Second we find values of the p line-shape parameters that are in agreement with PDG numbers. 
The mass is in agreement with the one obtained in e + e~ experiments. Furthermore the mass difference 
between the neutral and the charged ps is compatible with 0, so that no isospin violations are observed. 
The latter results improves the PDG values. 
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